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I describe a recent calculation (by me, Hoffmann, Liu, and Schaefer) of the chiral condensate in 
one-flavor QCD using numerical simulations with overlap fermions. The condensate is extracted 
by fitting the distribution of low lying eigenmodes of the Dirac operator in sectors of fixed topo- 
logical charge to the predictions of Random Matrix Theory. Our results are in excellent agreement 
with estimates from the orientifold large-N expansion. Much interesting physics surrounds this 
calculation, which I will highhght. 
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Theories like QCD, but with small differences, can teach us things about QCD. QCD with one 
flavor of dynamical fermion is such a theory: it is related, by a remarkable path, to the Nc—>-°° limit 
of ^ = 1 supersymmetric Yang-Mills theory, especially to one of its exactly-known observables, 
the gluino condensate. 

We lattice people are not the only ones interested in strongly coupled gauge field theories. 
Supersymmetry is a powerful tool for understanding these systems, and it is an important and active 
area of research, to extend results from supersymmetric theories to non-supersymmetric ones. One 
way to do this is to replace the degrees of freedom in supersymmetric field theories with new ones, 
while still preserving desirable features. The \aige-Nc (number of colors) limit is an important part 
of this program. (For an early attempt, see [||].) 

Recently, Armoni, Shifman, and Veneziano [^, ||, ^ ||] (ASV) suggested a new \aige-Nc ex- 
pansion with some remarkable features. In contrast to the 't Hooft laige-Nc limit ^ (Nc °°, g^Nc 
and Nf fixed, with quarks in the fundamental representation of SU {Nc)), quarks are placed in the 
two-index antisymmetric representation of SU {Nc)- Now in the Nc — > °°, g^Nc and Nf fixed limit of 
QCD, quark effects are not decoupled, because there are as many quark degrees of freedom as glu- 
onic ones, 0{NI) in either case. In Ref. [|5|] the authors have ai^gued that a bosonic sector of .JV = 1 
super- Yang-Mills (SYM) theory is equivalent to this theory in the large-A'c limit. = 1 SYM is a 
theory of adjoint gluons and their gluino (Majorana fermion) partners, and the equivalence of these 
theories in perturbation theory can be seen by comparing the vertices, as in Fig. |l], taken from Ref. 
[0]. The large-A^c QCD-like theory is called "orientifold QCD." 

The perturbative connection of orientifold QCD to ^ = 1 SYM is uncontroversial. In Ref. 
ASV have presented a nonperturbative proof of the connection. This proof has been extended by 
Patella ^ to lattice regularized theories. Recently, Yaffe and Unsal have argued that the proof 
of ASV is incomplete: that orientifold QCD and ^ = 1 SYM have a different phase structure 
on spacetimes with small compact dimensions, in which charge conjugation symmetry is sponta- 
neously broken. Only when the two theories have identical vacua can the proof hold. 

I am certainly not competent to comment more on this subject, so let us see what orientifold 
QCD might have to do with a lattice project: For Nc = 3, orientifold QCD is equivalent to QCD 
with a single quark flavor in the fundamental representation of SU (3). This equivalence can be seen 
in the first and second terms in the j8 function and in the lowest order anomalous dimension for the 
running quark mass (or quark condensate), as Table [T[ taken from Ref. [^] shows. This means that 
if the proof of nonperturbative equivalence is correct, nonperturbative quantities (in the bosonic 
sector) computed in super- Yang-Mills theory can be related to corresponding ones in one-flavor 
QCD, up to 1 /Nc effects. 

The analog of the quark condensate in ordinary QCD is the gluino condensate in .yV = 1 SYM. 
It can be calculated exactly in \wg&-Nc using saddle point methods (Ref. [|lO|] is a recent reference 
with a complete citation path). ASV used this exact result to estimate the quark condensate 
in one-flavor QCD from the value of the gluino condensate in SYM. They found (with our sign 
conventions) 

I = {0.014, 0.021, 0.028} GeV^ (1) 

in the MS scheme at jU = 2 GeV The spread of values gives their estimate of \/Nc corrections 
(basically l±l/3). 
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Figure 1: (a) The fermion propagator and the fermion-fermion-gluon vertex, (b) o/K = 1 SYM theory, (c) 
Orientifold daughter 

Now we come to the lattice: Last November, Schaefer and I had most of the solution to the 
problem of how to do Hybrid Monte Carlo for any Nf, using overlap fermions[pl|]. But what to 
do with it? I remembered the ASV prediction. Over Christmas vacation we started simulations to 
check this number and the preprint 1 12] came out in May. We hit their number bang on! 



We have to deal with a certain amount of imprecise language related to the condensate: The 
quantity (qq) is one definition of the condensate. Rather than measuring it directly, we will deter- 
mine the particular combination of the coefficients of the low energy effective field theory, £ = f^B, 
in the usual parameterization for Nf > 1 QCD, 

^2 = ^Trid^Ud^U"^) +B^Tr[M{U + U'^)]. (2) 

One expects that the quantity (qq) (as computed, for example, in a lattice simulation at some quark 
mass niq and simulation volume V) is a function of Z, /, rtiq, and simulation volume V. Z and / are 
the interesting quantities, and a direct lattice measurement of (qq) from several quark masses would 
have to be converted to a prediction of £ and /, by fitting it to the appropriate functional form from 
chiral perturbation theory. The same thing would have to be done if one measured observables, like 
the pseudoscalar mass and f^c, and used the GMOR relation to infer £. 

Nf = I QCD is a peculiar theory. Chiral symmetry is anomalous. There are no Goldstone 
bosons, just the eta-prime, which gets its mass through the anomaly. The £ which we are about 
to extract is therefore not an order parameter of spontaneous chiral symmetry breaking. How- 
ever, there still exists a well defined low-energy description of Nf = 1 QCD. It has been given by 



Leutwyler and Smilga [ ]13[ ] to which we refer the reader for details. They show that up to terms of 
order rrP'V the partition function is 

Z = exp [LVReime-'^)} (3) 

with 6 the vacuum angle. £ is the infinite volume zero quark mass limit of — (qq) at = 0. 

There is a reasonably inexpensive way to directly compute £: it involves measuring the distri- 
bution of the low eigenmodes of the Dirac operator, in sectors of fixed topology v in a simulation 
volume V. The probability distribution of individual eigenvalues A„ is given by Random Matrix 
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Table 1: A comparison of the beta functions and anomalous dimensions for the three theories 



Theory (RMT) 1 14, 15, 16] as a function of the dimensionless quantity X^LV, which depends para- 
metrically on the combination nigLV and, of course, Nf. We use the specific method and predictions 
from Refs. [jT^,|l8]l. 

The connection is only supposed to work deep in the e regime of QCD, but the eigenmode 
distribution is very robust and this does not seem to be a necessary constraint in practice. 
So the calculation involves several parts: 

• Perform simulations in sectors of fixed topology 

• Measure the eigenmodes and fit their distributions to the RMT expression: this gives the 
lattice-regulated a^H 



Determine a: we used the potential and the Sommer parameter [ |19| ] 



Convert to MS: we used the RI-MOM method [20] 



How we did the simulations Ql IQ is worth a paragraph. The last two items are pretty straight- 
forward. So is computing the eigenmodes. Fitting the distributions is also a little nonstandard, so 
I'll fill that in, too. 

When one is dealing with chiral symmetry on the lattice, it is very convenient to work with a 
fermion action which is chiral. That way, the physics of spontaneous symmetry breaking (including 
the anomaly) is not masked by explicit chiral symmetry breaking from the lattice discretization. 
While standard lore says that it is in principle possible to correct for chiral symmetry breaking 
effects in the bare action during the analysis, in practice this can be difficult and there is always the 
danger of new effects which one did not plan on (like exceptional configurations for Wilson-type 
actions). Why go looking for trouble? 



Thus we are led to use overlap [ ]21| , |22| ] fermions, which exactly encode chiral symmetry 
through the Ginsparg-Wilson [23] relation. It happens that, as an extra treat, it is possible to simu- 
late any number of flavors of overlap fermions, without requiring any degeneracy in the quark mass 
spectrum, using the exact Hybrid Monte Carlo algorithm. 

Due to the Ginsparg-Wilson relation, overlap fermions have the nice properties that the spec- 
trum of the squared massless Dirac operator = D^D commutes with 75 and has degenerate 
opposite-chirality eigenf unctions, apart from chiral zero modes and their partners at = 4Rq 
(where Rq is the radius of the Ginsparg-Wilson circle). The corresponding eigenmodes of H = 75D 



and D itself can be found by diagonalizing the 2x2 degenerate subspaces of H . Then [24, E3] the 
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contribution of the paired modes to the determinant of a single flavor of overlap fermions is given 
by detH^, the determinant of evaluated in a single chirality sector. It is included in the action 
with a single chiral pseudofermion for each flavor, AS = <paH(y{m)^^<p(y. The contribution of the 
zero modes can be included in the HMC by direct addition to the action at a topological boundary 
(the extra factor is Nf\Q\ ln(m^/ (2/?o)) for A^j degenerate flavors). For analysis, there is a different 
weighting of v = and v / configurations. 

What was never written down (maybe it is trivial) is how to initialize the pseudofermions: we 
need to generate random numbers to initialize The trick is to begin with a chiral random source 
R(y and to use the Zolotarev formula to construct (p^ = y/H(y{m)'^R(y. 

We used the algorithm of Ref. p6|]. In it, one must monitor the eigenmodes of the "kernel 
operator" h{—Ro) in D = /?()(1 + j5e{h{—R0))). When it develops a zero mode, the topology of the 
underlying gauge configuration changes and with it, there is a step discontinuity in the fermionic 
action. Random matrix theory wants eigenmodes in sectors of fixed topology. We generate those 
simply by forbidding tunneling events in the molecular dynamics evolution, and evolving in sectors 
of fixed V. 

To do the calculation, we collected data on a few P4's and P4E's for a few months, on 10^ and 
8^ X 12 lattices, at a lattice spacing of about 0.15 fm.. (The 10^ data set had about 500 trajectories 
each for winding number v = and 1.) 

Dynamical overlap is not cheap but it is completely feasible for small projects as long as 
one is willing to be creative (or maybe Baroque?) [27, 28]. The essential ingredient is a fat link 
gauge connection. We used stout links [29], three steps with p = 0.15. As far as the overlap goes, 
more smearing is better and 3 times 0.15 is about 2.5 times faster than the 2 times 0.15 of our 
previous work [27, ^]. The improvement comes from decoupling the fermions from UV gauge 
fluctuations which would generate small eigenmodes of h{—Ro). This decreases the condition 
number of h{—Ro) and speeds up the calculation of D. The dark side (if there is one) of a fat link 
action is that the fat links make the action more spread out than a conventional thin link action. 
Minimizing this spread is part of "action engineering," like minimizing the range of the fermionic 
couplings. We have performed the usual tests of locality on our action and never seen anything 
peculiar. Remember, thin links and fat links are just choices for the bare action which differ (in a 
Symanzik sense) through irrelevant operators. Formally, thin and fat link actions are both in the 
universality class of QCD. We are allowed to tune irrelevant operators as we please to ease the 
computational burden while preserving symmetries - not to do so is bad software engineering. 

If one wants to use eigenmodes, a kernel which looks "overlap-like" is also essential, otherwise 
the eigenmode part of the code is prohibitively expensive. (I have never been able to do anything 
with thin link Wilson-kernel overlap; it is too expensive for work stations. The problem is in the 
eigensolver: one begins with a set of trial eigenmodes which are then iteratively improved. Without 
some good idea, the beginning modes are typically random vectors. Improving them takes many 
iterations, which when done with the overlap action is very slow. One can gain a lot of time using 
better eigenmodes, from some "overlap-like" action, as an intermediate step. The overlap actions 
I use are built on "overlap-like" kernels, but the overlap action with a Wilson kernel has a very 
different spectrum from a Wilson action.) 

All of this is well documented for quenched simulations [pO[]. I can't help thinking that there 
are more tricks out there. Oh, for another factor of five speedup... 
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So we collected a set of fermionic eigenmodes. We want to fit their probability distributions 
to RMT formulas. The analysis is a little nonstandard, because the eigenmodes are continuously 
distributed. A referee led us to the Kolmogorov-Smirnov test [ pi] ] as a measure for the goodness of 
the fit. It compares the cumulative distribution function of the data C(x) to the theoretical prediction 
P{x) = J^^f{x)dx. C{x) is the fraction of eigenvalues with a value smaller than x. 

The quantity of interest is the largest deviation of P and C. D = maxj^ |^(-^) —C{x)\. From this 
the confidence level is given by 



Qks [{Vn + 0.12 + 0.n/VN)D) (4) 

with 



e^s(A) = 2X(-)^--iexp(-2/A2). (5) 

;=i 

In fits to a single eigenmode distribution we maximize this quantity. When fitting to more than 
one mode, we maximize the product over the individual confidence levels. The errors on the fit 
parameter £ are determined by the bootstrap procedure. An example of such a fit is shown in 
Fig. ^. After a lot of angst (which modes to fit, what about correlations...) Za^ turned out to be 
remarkably robust: it didn't matter what we did. 

Completing the calculation with the lattice spacing from the Sommer parameter and the match- 
ing factor from RI-MOM, we found 

rlZ(MS,n = 2 GeV) = 2,(^1, a) x la^ x {'^f 

= (0.86(3)) X (0.0096(3)) x (3.37(10))^ 

= 0.317(22). (6) 

With ro = 0.5 fm, this is 

L{MS, PL = 2 GeV) = 0.0194(20)GeV3 (7) 

which agrees pretty nicely with Eq. [T[ 

The summary of McNeile [ ]32| ] shows that the Nc = 3 condensate is not very Nf dependent. 
Indeed, Schaefer, Liu and I just finished [33] an Nf = 2 measurement using basically identical 



techniques to what I have described here, and we find I = 0.0225(25) GeV^ =(282(10) MeV)^ 
There is actually an annoying systematic in this number: In finite volume, there is a first order 
correction to the condensate, basically the one loop graph from Goldstone bosons which are emitted 
from the propagating pseudoscalar and which are absorbed at an image point of the vertex. The 
correction is £ — > pZ, where 

and c{li/l) depends on the geometry[34]. This term is absent in Nf = 1 because there are no 
Goldstones. We haven't measured /;;:, but with 93 MeV, in our volume for Nf = 2, p ^ 1.43. 
Fortunately, people publish not £! 

After the paper came out, Veneziano reminded me that their prediction was really for Z/A^ 
"without going through actual experimental numbers." For us lattice people, a ratio of nonpertur- 
bative quantities with no (or minimal) intrusion of perturbation theory is much cleaner than T/A^, 
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Figure 2: Cumulative distributions from RMT fit to eigenmode distributions. For each mode, the three 
dotted lines, running from left to right, correspond to fits to the lowest mode in the v = 1 sector only, the 
lowest mode v = and |v| = 1, and the lowest mode in v = only. 



and the definition of A is exquisitely sensitive to a determination of a coupling: recall that 



Tadpole-ology (plaquette —>■ ay {q* ) ^ A) at one value of the lattice spacing is just too unstable 
to be useful. However, the ALPHA collaboration has used the Schrodinger functional real space 
renormahzation group to compute [3^, 36] the quantities 



A{MS)ro = 0.62(4) (4); Nf = 2 
= 0.60(8). A^/ = 



(10) 



If we take A{MS)ro = 0.61(6) for Nf = 1, then 



= 1.4(4) (11) 

while ASV want 0.6 to 1.1 for the ratio. This is not so bad! (The difference between this and Eq. |] 
is a 15 per cent shift between their A and the interpolated ALPHA value.) 



Armoni, Shore, and Veneziano have also predicted £ for other N/s p7[]. They add Nf — I 
fundamental flavors to the mix, and get 

for the RGI condensate. They need a coupling constant to convert this to an MS number. Their 
publication only presents a band, since when they do the conversion from RGI to MS they consider 
a range of coupling constants. However, taking the Nf = 2 value of K{MS) from Eq. |T^ and 
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inverting it to give a coupling constant, Z'/^ = 247 MeV. More accurate lattice measurements of £ 



(and A) vs Nf could test the Nf dependence in Eq. [12. 

I am not sure what to do next. Unfortunately, the only nonperturbative quantity which can 
be computed in SYM is the gluino condensate (as far as I can tell, from asking many people), 
although to compensate, it is supposedly exact. So other predictions typically involve ratios of 
masses without some connection to the condensate: 

2 

mt., 

-JL = i + 0{l/N,) (13) 



(From Ref. [g]; the 0{\/Nc) corrections have been computed by Sannino and Shifman [)38[]), 
and a similar degeneracy for hybrids [^9|]. An Nf = \ meson is just like a flavor singlet meson 
in ordinary multiflavor QCD, so both the r]' and a have disconnected ("hairpin") contributions. 
These are difficult and noisy. In addition, the scalar operator has a vacuum expectation value, so 
its signal is like a scalar glueball's: the exponential which gives the mass dives under a constant 
background. I have successfully avoided trying to do these for four months, now. 

Nf = 1 QCD has other intriguing properties: At negative quark mass, Nf = 1 QCD may have a 
phase in which CP is broken. (This can occur for any Nf if the quark mass matrix has positive and 
negative eigenvalues. This observation goes back to Dashen [pO|]; see Smilga [41] and Creutz [42] 
for recent discussions.) It is possible to simulate QCD at mq < (or even complex mass) with 
overlap fermions by reweighting a real-mass simulation. The relevant derivation has been given by 
Diirr and Hoelbling [p3|], who briefly studied the Schwinger model. Complex mass is equivalent to 
a theta vacuum, another unvisited area of QCD for lattice simulators. 

"The purpose of computing is insight, not numbers" (Hamming), so did we learn anything? It 
is always easiest to say No. The only quantitative prediction from ^ = 1 SYM is the condensate, 
and it has unknown 1 /Nc corrections. The spectrum of QCD depends on Nf in an uncontrolled 
way. The only way to make predictions relevant to the real world is with simulations with three 
flavors of light quarks, all in the chiral regime. Even then, one must simulate at physical quark 
masses, unless the observable being studied has a well-behaved expansion in chiral perturbation 
theory allowing one to extrapolate in quark masses. 
And yet- 
When we teach about the spectrum of hydrogen in an introductory quantum mechanics class, 
the story is not "we do the calculation and the answer is 13.6 eV." There is systematic expansion 
(in a or equivalently in v/c) which allow us to make successively more accurate predictions. The 
parameters we have at our disposal in QCD are Nc, Nf, the color representations of the quarks, and 
the quark masses. Perhaps the zeroth order QCD calculation (like V{r) = —e^/r for hydrogen) is 
some extreme value of one or all of these parameters. Lattice tests of QCD-like theories might 
tell us new things about QCD, if they could validate extrapolations of analytic results from those 
theories. We won't know if we don't try! 

And this test worked: ASV successfully predicted the Nf = I condensate. 
I would like to thank my collaborators Roland Hoffmann, Zhaofeng Liu, and Stefan Schaefer 
for many conversations, and I am grateful to Adi Armoni, Francesco Knechtli, Francesco Sannino, 
Misha Shifman, Matt Strassler, Mithat Unsal, Gabriele Veneziano, and Larry Yaffe for correspon- 
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dence. And thanks to the organizers for putting on a fantastic conference! This work was supported 
by the US Department of Energy. 
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